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ABSTRACT 

A Newtonian model of non-conductive, charged, perfect fluid tori orbiting in combined spherical gravita- 
tional and dipolar magnetic fields is presented and stationary, axisymmetric toroidal structures are analyzed. 
Matter in such tori exhibits a purely circulatory motion and the resulting convection carries charges into per- 
manent rotation around the symmetry axis. As a main result, we demonstrate the possible existence of off- 
equatorial charged tori and equatorial tori with cusps enabling outflows of matter from the torus also in the 
Newtonian regime. These phenomena represent qualitatively a new consequence of the interplay between 
gravity and electromagnetism. From an astrophysical point of view, our investigation can provide insight into 
processes that determine the vertical structure of dusty tori surrounding accretion disks. 
Subject headings: accretion, accretion disks — magnetic fields 



1. INTRODUCTION 

As demonstrated by a variety of observational phenomena 
ranging from spectral and timing properties of accreted gas 
in different energy bands, luminosity profiles, and velocity 
dispersions of stars forming the nuclear clusters, up to direct 
measurements of stella r motion in the Galactic center ( |Kro-| 
[Ek 2004; Eckart et al.]|2005), supermassive black holes are 
frequently presenPm the nuclei of galaxies. The radio mor- 
phology of jets emerging from galactic nuclei and the skewed 
double-horn profiles of many X-ray and UV spectral lines 
suggest that these objects obey axial rather than spherical 
symmetry (Urry & Padovani 1995). As a working hypothe- 
sis, a central black hole with a radiating inner accretion disk 
is thought to be surrounded further out by a torus of obscuring 
material (dust). 

Typical masses of supermassive black holes fall within the 
range of M. ~ 10''-10^Mq. The central mass is related to 
the gravitational radius of the black hole, 7?g = GM,/c^ ~ 
1.5 X W^(M,/\(fMQ) cm ~ 10"^ pc. On tlie other hand, 
the typical radius of the torus, as derived from observations, 
greatly exceeds the gravitational radius, ~ 10-100 pc. The 
mass of the torus is orders of magnitude smaller than the black 
hole mass, Md <C M,. Therefore, is usually neglected 
in discussions, however, it may be important for the internal 
structure of the torus. In fact, the dusty tori are thought to be 
on the verge of self-gravitational instability, which can lead 
to the fragmentation of the torus when T oomre's criterion of 
stability is violated (Cofli n & Zahn|2008[ ). 

Different types of active galactic nuclei in Seyfert galaxies 
can be unified by introducing some form of obscuring tori, 
which are b elieved to encircle the central blac k hole ( |An- | 
tonucci & M iller 19^ |Urry & PadovMiil[l995] l. The pres- 



ence of a geometrically and optically thick dusty st ructure is 
an essential component of the unification scheme (Honig & 



motion and its rotation defines the symmetry axis of the sys- 
tem, which may or may not be aligned with the rotation axis 
of the black hole. It is generally expected that the alignment is 
enforced at least in the inner regions by the Bardeen-Petterson 
effect (Kumar & Pringle 1985). 

Equilibriu m figures of gaseous tori have been studied in 
great detail ( Kozlowski et al.|1978[ Abramowicz et al.|1978 



Kato et al. 20081, however, the vertical component of the 
pressure gradient, required to maintain the equilibriu m, does 
not seem to be sufficient in dusty tori (see, e.g.. Mur phy &] 
,Yaqoob (2009) and references cited therein). Therefore, it ap- 
pears that these structures have a tendency to collapse toward 
the plane. Moreover, the self-gravity influences the structure 
of massive tori (as well as geometrically thin Keplerian disks) 
once the local mass density of the ac creted gas at a given lo- 
cation exceeds ^ ri^/47rG, see, e.g., Shlosman & Begehnan 



|Kishimoto 20 lO) . The torus structure is thought to be inho- 
mogeneous, in the form of molecular/dusty clumps (clouds). 
The matter forming the torus follows a common bulk orbital 



( |1987| ); Hure ( 1998). Despite the fact that signahires of obscu 
ration (especially those seen in X-ray spectra) and variability 
properties strongly indicate the need for a significant vertical 
extent of obscuring tori in many Seyfert type 2 galaxies, the 
physical model for the tori remains uncertain. 

The above-given arguments suggest that the vertical struc- 
ture of dusty tori needs fur ther discussion. For examp le, it 
has been proposed recently ( [Czerny & Hryniewicz|201 1) that 
vertical motions of the dust clumps play an important role. 
These can be driven by radiation pressure in the region where 
dust can survive evaporation. In fact, the irradiation in active 
galactic nuclei is very intense and it presumably leads to the 
dust evaporation below a critical radius (typically at a fraction 
of parsec) where the dust temperature reaches about lO-'K. 

It is important to note that the irradiation of dust particles 
should produce a posit ive net electrical charge by photoion- 
ization (Ho ranyi 1996 Vladimirov et al.|200 5). On the other 
hand, plasma electron and ion currents are continuously enter- 
ing the grain surface, so the sign and magnitude of the equi- 
librium charge depend on the total currents that are absorbed 
and emitted from the grain surface. This is a complex process 
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Figure 1. Region of pai'ameters K2 and rc , enabling the existence of fam- 
ily I and II equatorial tori with a uniform distribution of the specific angular 
momentum, l(r, 9) = K2 = const. The dashed curve divides the region into 
subregions with positively/negatively (+/— ) charged tori. The interval of pos- 
sible values of rc for the choice A'2 / V GM = 2.5 is marked and enlarged. 



that depends on many parameters; simultaneously, it appears 
to be an important mechanism that can govern the torus struc- 
ture. The electrostatic charge is one of the essential parame- 
ters that controls the dynamic s of dust gra ins embedded in the 
surrounding cosmic plasma (I shihara|2007| l. 

When the electromagnetic force, due to the pervasive ex- 
ternal magnetic field, is taken into account, an electrically 
charged dust can establish vertically extended structures that 
'levitate' above and under the equatorial plane. Namely, 
charged test particles can be bounded within off-equatorial 
wells of effective potential, which suggests the existence 
of very dilute halo toroidal formations consisting of non- 
interacting, charged grains, ions, or electrons. The Newto- 
nian study of charged dust grains orbiting in planetary mag- 
netospheres an d forming halo orb its was published, e.g., in 
Howard et al.| (|l999 ); Dul Un et aL] 12002), while the question 
of whether sucn halo orbits can appear also in strong gravita- 
tional field s near compact objects was positively answered in 
Kovar et al. (2008); StuchliTcet al. (2009); Kovaret al. (20lO]l. 



Of course, in many astrophysical scenarios such simple test- 
particle approaches fail because of higher densities of charged 
matter in reality. Then, possible approaches follow from the 
kinetic theory (suitable for lower density matter) or from the 
hydrodynamics (suitable for higher densities). 

In this paper, we apply the Newtonian hydrodynamical ap- 
proach, modeling perfect fluid tori with electric charge spread 
through the fluid of infinite resistivity. This model of dielec- 
tric fluid represents an opposite limit to the well-known ideal 
magnetohydrodynamics with zero resistivity as a goo d ap- 
proximatio n of many astrophysical plasmas, see, e.g., |Pun-| 
|sly| ( |2001[ ). The general relativistic description of charged 
perfect-fl uid disks of infinite resistivity was published re- 
cently in [Kovaf et al. (2011 1, where the tori encircling static 
electrically charged Reissner-Nordstr0m black hole were an- 
alyzed. Note, however, that the introduced Newtonian model 



is basically (technically) simpler and more convenient for our 
primary investigations of halo tori formed far away from the 
central black hole (i.e., in a weak gravitational field) due to 
the presence of external magnetic fields. 

We approach the problem within the framework of a two- 
dimensional (axially symmetric) toy model of electrically 
charged tori rotating in spherical gravitational and dipolar 
magnetic fields. By constructing the equilibrium figures, we 
find the constraints on the spatial distribution of charge and 
matter densities, and on the corresponding angular momen- 
tum profile. The main aim of our investigation is to explore 
how the well-known structures of electrically neutral gaseous 
tori become modified by the interaction of (small) charge den- 
sity with the poloidal magnetic field, leading to vertically ex- 
tended steady distribution which can reach further out of the 
equatorial plane and remain in a stable permanent azimuthal 
circulation. 

The paper is organized as follows. The equation of motion 
for a perfect fluid with spatially distributed electric charge 
(Euler equation) and the corresponding condition of hydro- 
static equilibrium are formulated in Sec. |2] In Sec. [3] this 
condition is analyzed for incompressible nuid and possible 
stationary toroidal configurations are found. The relevance 
of the results obtained for a polytropic fluid is discussed in 
Sec.|4] InSec.|5] the use of Newtonian formalism for compact 
objects is demonstrated by replacing the Newtonian gravita- 
tional potential with the Paczyriski-Wiita one and the corre- 
sponding modifications of isobars are discussed. Concluding 
remarks are given in Sec.|6] 

2. NEWTONIAN MODEL 

The equation of motion for a perfect fluid (no viscosity) or- 
biting in g ravitational and electrom agnetic fields, i.e., Euler's 
equation ( |Landau & Lifshitz|I987 1, has the form 



QmidtVi + vj^ jVi) = -^iP- gm^i^ + gei.Ei + eijkvjB''), (1) 

where and are the mass-density and charge-density, re- 
spectively, P denotes the pressure, v is the velocity field in 
the fluid, and $ corresponds to the gravitational potential. An 
electromagnetic field is described by its electric part E and 
magnetic part B. The last term on the rhs of Euler's equation 
(fill thus corresponds to the Lorentz force density. 

In our analysis, we assume stationary and axisymmetric 
flow of test charged perfect fluid in external spherical grav- 
itational and dipolar magnetic fields. In spherical polar coor- 
dinates (r, 9, (p), 



GM 



£,=0, i = (r, 
cos 9 



sm9 

Br = 2lJ, r— , Bg = ^ — r— , 



(2) 
(3) 
(4) 



where M is the mass of central object, i{r, 9) is the specific 
angular momentum of fluid elements, and /i > corresponds 
to the magnetic dipole moment of the external magnetic field. 

Conservations of mass and electric charge are described by 
corresponding continuity equations 



9,£'m + V,(£lmV') = 0, 
5,& + V,(£leV') = 0. 

For stationary and axisymmetric flow. 



(5) 
(6) 

(7) 
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Figure 2. Positively ciiarged equatorial torus (family /). The meridional section and equatorial profile of the pressure field (left) and the specific electric charge 
(right) are shown. 



thus, both the continuity equations (jSj) and (|6]) are fulfilled 
automatically. 

The proposed fluid dynamics is described by two partial 
differential equations, following from Euler's equation ([T]i: 



use the limit of incompressible fluid, for which g,^ = const. 
Equations ([8]l and (|9]l then have the form 



GM 



dP _ 
dr 

\dP ^ 
r ov r 



smt 



COSf 



(8) 



(9) 



1 dP 
Qm dr 
1 dP 

Qm 99 



GM V 



V 



sm 



Hqv^— \ =A{r,e), (11) 



B{r,9). (12) 



In order to solve this set of equations, it is useful to assume 
charge density in the form 



Qe = gmq(r, 0), 



A necessary and sufficient condition for the existence of the 
solution of the above-mentioned set of partial differential 
equations, describing a distribution of pressure in the fluid, 
is given by the integrability condition 



(10) 



where the function q{r^ 9) describes the specific charge distri- 
bution in the fluid. To complete the set of equations, we need 
to establish an equation of state. 

3. INCOMPRESSIBLE FLUID 

When there is no significant compression or expansion of 
the fluid (liquid or gas) anywhere along flow Unes, we can 



dA_dB 



(13) 



having the form 
Idv 



2v„ 



r d9 



^ cot6' ^ 



dr 



3qv^ 



cost 



sinfc 



dq dv^\ ^cos9 ( dq 9v„' 



'89 ^ 89 



dr dr 



= 0.(14) 
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Figure 3. Negatively charged equatorial torus (family I). The meridional section and equatorial profile of the pressure field (left) and the specific electric charge 
(right) are shown. 

The integrability condition has the form of the partial differ- 
ential equation for two unknown functions y^(r, 6) and q{r, 9). 
In order to find their proper forms, it is necessary to spec- 
ify somehow at least one of these functions. In this paper, 



general forms of such distributions satisfying the integrabiUty 
condition ( [T4] l: 



we choose the function v,^(r, 9), i.e., the way of rotation of 
the fluid. Analysis of condition ([14]) reveals that in the case 
of q = (oT II = 0) the orbital velocity v^p can be written in 
a separate form v^(r, 9) = v(r)u(0), where the functions v(r) 
and u{0) satisfy differential equations (ensuring that the first 
parenthesis in the integrability condition ( [T4| is zero) 



g(r,0) = Cr-3(^'-'>/2, 
q{r,9) = Cr^l^(jim9y^'^' 
q{r,9) = Cr-^'^'l^ sin^6', 
^(r, 6') = C(sin6')^*'-^''; 



(17) 
(18) 
(19) 
(20) 



dlnv dlnw 
dlnr dln(sin6') 

with Ki being a separation constant. Solving Eqs. ( [T5] l 
obtain a general form of the orbital velocity profile 

v^{r,e) = K2{rs,m9f\ 

where K2 is an integration constant. We will use such general 
form of v^p{r^ 6) also in situations, in which qy^Q (and fi y^Q). 

Further, we assume that also the distribution of a specific 
charge in the fluid, i.e., the function q(r, 9), can be given in 
the separable form q(r, 9) = qiir)q2(9)- Then there are four 



(15) 



we 



(16) 



C is a constant of integration. Further, we shall refer to distri- 
butions ([T7]i-([20li as families I, II, III and IV, respectively. We 
can see that for certain values of constant Ki, there are some 
overlaps between those families: for Ki = 1, famiUes I and IV 
coincide, for Ki = 0, families I, II as well as III, IV coincide, 
and for Ki =-l this is true for families II and III. 

In the center of the torus {rc,9c), the pressure takes its 
maximal value. The necessary condition VP = (especially 
dP/d9 = 0, see relation ( 12 1) points out that the pressure max- 
ima can be located either tn the equatorial plane {9^ = 7r/2) or 
even out of the equatorial plane (9^^-k/2), where the possible 
off-equatorial maxima are supported by the specific charge of 
the torus with its distribution along the off-equatorial circle 
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Figure 4. Positively cliarged equatorial torus (family II). The meridional section and equatorial profile of the pressure field (left) and the specific electric charge 
(right) are shown. 



given by the relation 



From Eq. ( 11 1, it follows that the orbital velocity v^{r, 9) > 
is given by tne relation 



(21) 



The condition VP = is only a necessary condition for an 
existence of local extrema. To verify thek real existence, we 
must evaluate the second order derivatives and construct the 
Hessian matrix 



IGM 



AGMr^J 2rl ^ ' 



there. 



Analysis of conditions ( 23 1 shows that equatorial local 



d-P d^P 

-TJ - I df- drdB 

\ d-P ifp 

dear 



(22) 



maxima of the pressure exist for all four families (17i-(20i 
of the specific charge q{r, 6) if, in the center of torus, tne con- 
stants K\,K2, and C satisfy the conditions 



in the loci of pressure maxima, two conditions must be satis- 
fied: 

-gj \r=r,.e=e^ < 0, detn\,.=r^.e=e, > 0. (23) 



3.1. Equa to ria I to ri 

A local maximum of the pressure corresponding to the cen- 
ter of the torus is located in the equatorial plane, 9c = ir/l. 



family I: 



K 



2_ J/Ti+l < ■'^l 1 



CM 
GM'" 



3= 



for/Ti ^ — , 
5Ki + l < 5' 



^ 3' 



/iC = K2rl 



,(5^,+l)/2_ 



GM 



(25) 

(26) 
(27) 
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Figure 5. Region of parameters K2 and , enabling the existence of fam- 
ily III equatorial tori with the specific angular momentum profile £(r, 8) = 
/f2(rsine)'/3. The dashed curve divides the region into subregions with pos- 
itively/negatively (+/-) charged tori. The interval of possible values of for 
the choice K2/ V GM = 2.0 is marked and enlarged. 

family II: 



GM " 3' 



^lC = K2r\ 



,(2/ri+i)/2_ 



GM 

1 



-(2/ri+i)/2 



family III: 



2a:,+i . _ 
GM ' 3 



Ki+2 
5Ki+4 
2 



,(5/fi+4)/2_ 



GM 



.(Ki+2) 12 



family IV: 



Jf2 

■^2 p.Kx + \ 

GM ' 



/r2 

■^2^2A:i + 1 ^ 



Ki+2 
2 



for Ki ^ -2, 



GM'" ' 3' 



GM 

1 

K2 



(28) 
(29) 
(30) 

(31) 
(32) 
(33) 

(34) 

(35) 
(36) 



In the following paragraphs, we present illustrative exam- 
ples of equatorial tori for each of the specific charge distri- 
bution families. A particular solution of the set of equations 
( [TT| and ( [T2I ) is obtained when the constants Ki, K2, and C are 
specified with respect to conditions ( 25 1-( 36 1 for a given fam- 
ily. Instead of the value of C, the location of the torus center 
can be prescribed. Of course, the choice of constants should 



be physically meaningful, reflecting realistic ideas about the 
motion of matter Such a basic idea is, e.g., subluminal mo- 
tion inside the torus, v^p < c (c is the vacuum speed of light), 
which is not automatically incorporated in the Newtonian de- 
scription. 

3.1.1. Family I equatorial tori 

As an example, we choose Ki = -1, which corresponds to 
the torus with uniform distribution of the specific angular mo- 
mentum, £{r, 6) = K2 = const. In this case, the specific charge 
distribution and the pressure field are given by the relations 



q{r,e) = Cr\ 



P(r,0) = PQ+g, 



where 



GM 



gmvl, + gmK2tiCln 



K2 



rsinO 



(37) 



(38) 



(39) 



and Pq = const is an integration constant, determining the sur- 
face of zero pressure (boundary of the torus) and, thus, the 
torus extension and thickness as well. 

To obtain a particular pressure distribution, we further 
choose Pq = 0, Ko/VGM = 2.5, and, in accordance with 
conditions p5]l a nd (26i graphically presented for the case 
A"! = -1 in Fig.yj = a Further from relation ( [27| we obtain 
[iCjyJ GM = 2.78 X 10"^, which corresponds to the positively 
charged torus (C > 0). The meridional section and equato- 
rial profile of the resulting pressure field, as well as the spe- 
cific charge distribution, are presented in Fig. 2 ' We see that 
there are toroidal isobaric surfaces corresponding to equilib- 
rium stationary fluid torus with the spherical distribution of a 
positive specific charge. 

The negatively charged torus is determined by the condition 
C < (C is given by relation ( 27 1) which, together with con- 
ditions ( 25 1 and ( 26 1, puts a limit on the location of the torus 
center. For the given values of Ki = -1 and K2/V GM = 2.5, 
the center of the negatively charged torus must be located in 
between the radii 6.25 and 9.375, see Fig.[T| We choose = 8, 
which gives iiC/y/GM = -0.012. In this case, the equatorial 
pressure profile has two local extrema (maximum and min- 
imum) around which the toroidal isobaric surfaces are cen- 
tered. Nevertheless, only those around local maximum corre- 
spond to toroidal fluid body with pressure gradient decreasing 
from the center to an outer edge, which is, in this case, given 
by the critical isobaric surface. Pa-it = 0, with two cusps out 
of the equatorial plane (Pq = -0.13 in this case). The cusps 
enable outflow of matter from the torus to an outer space 
due to the violation of hydrostatic equilibrium (the so-called 
Paczyhski mechanism known in the theory of accretion disks). 
The meridional section and equatorial profile of the pressure 
field, as well as the specific charge distribution, are presented 
in Fig. [3] 

3.1.2. Family II equatorial tori 

Again, we choose Ki = -1. The corresponding specific 
charge distribution and pressure field are described by the for- 
mulae 

(40) 



sin 



In all figures, units G = M = I are used. 
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Figure 6. Positively charged equatorial torus (family HI). The meridional section and equatorial profile of the pressure field (left) and the specific electric charge 
(right) are shown. 

3/2 



GM I 2 2 

P(r, 9) = Po+ 1^0mV^+- QmKl^lC 

r 2 ^3 



where 



K2 

rsin^' 



(41) 
(42) 



K2 = £ = const, and Pq = const. 
Further, we choose Pq = -0.02, K2/V GM = 2.5 and, in ac- 



(maximum) and the pressure field has the same topology as in 
the positively charged case. 

3.1.3. Family III equatorial tori 

Here, we choose Ki = -2/3^ The corresponding specific 
charge distribution and pressurerield are described by the for- 
mulae 

(43) 



q(r, 9) = Cr sin' 9 



cordance with condition (29 1 depicted for the case Ki = -1 in 
Fig.[l] rc = 6. According to relation ^C/VGM = 0.041. 
The meridional section and equatoriaTprofile of the pressure 
field, as well as the specific charge distribution, are presented 
in Fig.|4] The relevant part of the pressure field has a toroidal 
topology, while the specific charge distribution has a cylindri- 
cal topology. The presented case corresponds to a positively 
charged torus. 

The negatively charged torus with the same values of Ki 
and K2 would be centered around a circle with the radius 
6.25 < Tc < 9.375, see Fig. [T] The equatorial pressure pro- 
file of negatively charged ton still has only one local extreme 



P(r,9) = Po + g, 



where Pq = const. 



GM 3 



5/3 



7gmV^ + -gmK2lJ.C 



sm 



v„ = 



K2 



(r sin 61)2/3' 



(44) 



(45) 



^ The choice Ki = —1 used in previous examples gives the specific chai'ge 
distribution of t he last example and thus leads to the pressure field described 



by formula j41| . 
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Figure 7. Negatively charged equatorial torus (family III). The meridional section and equatorial profile of the pressure field (left) and the specific electric charge 
(right) are shown. 

and K2 = const. The distribution of the specific angular mo- 
mentum is given by the relation 



1/3 



(46) 



Further, we choose K^jsJ GM = 2.0 and, in accordance with 
conditions ( [3T| ) and ( [32| i graphically presented for the case 
Ki = -2/3 inFig. |5] rc = 12. These choices give positively 
charged torus, since according to relation (33i ^C/^/GM = 
1 .958. The structure of isobars, the equatoriaTradial profile of 
the pressure, and the specific charge distribution in the fluid 
are shown in Fig.|6] We see that there is one self-crossing iso- 
bar of the 'Roche^obe' type surrounding closed isobars and 
forming a cusp that enables the outflow of fluid from the torus 
onto the central body. Note that the constant Pq is chosen to 
obtain P^it = for this critical isobar The location of the cusp 
corresponds to the position of a local minimum in the radial 
pressure profile. The surfaces of constant specific charge have 
cylindrical topology. 

In order to obtain a negatively charged torus with the same 
values of Ki and K2 as in the case of positively charged torus 
presented above, we must shift a torus center to the position 
64 < < 216, see Fig. |5] We choose rc = 80, which gives 



fiC/VGM = -0.665, and Po = -0.0015. The resulting pressure 
and specific charge distributions are shown in Fig.|7] There 
is no self-crossing isobar, only the toroidal ones exist, corre- 
sponding to a negatively charged toroidal structure with the 
cylindrical distribution of the electric charge. 

3. 1 .4. Family IV equatorial tori 

As an example, again we choose tori with uniform dis- 
tribution of the specific angular momentum, £{r, 9) = const, 
which are characterized hy K\ = -1. The corresponding spe- 
cific charge distribution and the pressure field are described 
by the formulae 



^(r, 61) = C sin" 61, 



P(r,0) = Po + ft 



where 



GM 



rsmt 



(47) 



(48) 



(49) 



K2 = i = const, and Pq = const. 
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Figure 8. Region of parameters K2 and r^, enabling tlie existence of family 
IV equatorial tori with a uniform distribution of the specific angular momen- 
tum, £(r, 8) = K2= const. The dashed curve divides the region into subregions 
with positively/negatively (+/— ) charged tori. The interval of possible values 
of Tc for the choice & / V GM = 2.5 is marked and enlarged. 



pressure distribution, we choose 
rc = 4 in accordance with condi- 
(35 1 represented for the case Ki = -1 in 



For a particular 
K2/VGM = 2.5 ^d 
tions ( [34[ ) and 

Fig. [8] These choices give a positively charged torus, since, 
according to relation (36i, iiC/\/GM = 3.6. The structure of 
isobars, the equatorialradial profile of the pressure, and the 
specific charge distribution in the fluid are shown in Fig. |9j 
As in the case of family III tori, there is one self-crossing 
isobar of the 'Roche-lobe' type surrounding closed isobars 
and forming a cusp (the constant Pq is chosen to obtain 
^ciit = for this critical isobar). The location of the cusp 
is given by the position of the local minimum in the radial 
pressure profile. The surfaces of constant specific charge 
correspond to vertical planes. 

The negatively charged torus with the same values of Ki 
and K2 as the positively charged one presented above must be 
centered around a circle with the radius 6.25 < < 9.375, 
see Fig. [s] We choose = 8, which gives /iC/ %/ GM = -5.6, 
and Pq = -0.02. Resulting pressure and specific charge dis- 
tributions are shown in Fig.[TO] There is no self-crossing iso- 
bar, only the toroidal ones exist, corresponding to a negatively 
charged toroidal structure with the radial distribution of the 
electric charge. 

3.2. Off-equatorial tori 

Due to the mutual interaction between electric charges and 
the magnetic field, the torus center can be lifted away from 
the equatorial plane, i.e., 9^ ^ 7r/2. The precise location of the 
center is governed by the interplay with pressure gradients. 
Since the specific charge is given by relation pT] ) there, we 
immediately see that for positively/negatively rotating torus 
(v^ > Q/v^ < 0) and the positive/negative direction of ex- 
ternal magnetic field (/i > 0//i < 0), the off-equatorial torus 
must possess a negative electric charge. On the other hand, 
positively/negatively rotating tori in negatively/positively ori- 



ented magnetic fields possess a positive electric charge. Next 
we will focus only on off -equatorial tori with V;^ > orbiting 
in the external magnetic field characterized by /i > 0, i.e., on 
the negatively charged off-equatorial tori. 

In the torus center, there is VP = 0. Adding relation (21 
for the specific charge to Eq. (fTTli and putting dP/dr = 0, we 
obtain the relation for the orbitaTvelocity in the center 



v^(rc, i 



2GM 



(50) 



Consequently, relation (21 1 for the specific charge in the cen- 
ter takes the form 



1 

2/i 



2GM 



3r„ sin^ 



(51) 



Analysis of conditions ( 23 1 shows that off-equatorial local 
maxima of the pressure can exist for families II, III and IV of 
the specific charge distribution q(r^ 6). In detail, in the torus 
center the constants Ki, K2, and C must satisfy conditions 



family II: 



/jlC 



1 

"2' 
K2 



i2Ki 



4K,-1 



family III: 



/Ti €(-3.26,-0.41), 

^C = -^r^^'''^'>'\sme,f'-'. 



famOy TV: 



K, <- 



10-2 



-0.387, 



K2 



.(Ki+2) 



(sin6'c)' 



MKi-l) 



(52) 
(53) 

(54) 
(55) 

(56) 
(57) 



Further, for each of these three families, we present an illus- 
trative example of off-equatorial toroidal structures. Again, 
we must specify the integration constants Ki , K2, and C. How- 
ever, instead of C, it is more convenient to choose the latitude 
of the torus center, i.e., 6^. 

3.2.1. Family II off-equatorial tori 

In order to be consistent with previously studied equato- 
rial tori, we choose the family II tori with uniform distribu- 
tion of the specific angular momentum, £{r, 9) = K2 = const. 
In this case, Ki = -1 and the specific charge, pressure and 
orbital velocity fields are given by relations ( 40 1-( 42 1. Fur- 
ther, we choose Pq = -0.015, K2/y 



GM = 2.5, and6'c = ±7r/4 
(as the pressure field is symmetric with respect to the equa- 
torial plane). The radial coordinate of tori centers is deter- 
mined by relations ( [42] i and ( [SO] ) for orbital velocity, which 
have to be fulfilled simultaneously in the center The cho- 
sen values of K2 and 9^ give r^ = 18.75 and, according to 
GM = 



(53 I, fiC/VGM = -1.633. The meridional sections and off- 
equatorial (9 = 7r/4) profiles of both the pressure field and 
specific charge distribution are presented in Fig.fTT] showing 
toroidal off-equatorial structures above and underthe equato- 
rial plane with cyUndrical distribution of the specific charge. 
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Figure 9. Positively cliarged equatorial torus (family IV). The meridional section and equatorial profile of the pressure field (left) and the specific electric charge 
(right) are shown. 



3.2.2. Family III off-equatorial tori 

As in the case of equatorial family III tori, we choose 
Ki = -2/3. The specific charge distribution, pressure field, 
orbital velocity field and specific angular momentum distribu- 
tion are thus given by relations ([43]l-([46]l. Further, we choose 
K2/V GM = 2.0 and 9^ = ±7r/4, giving the radial coordinate 
of tori centers = 864. Subsequently, [iCjyjGM = -48. For 
Pq = -2.5 X lO"'*, the meridional sections and off-equatorial 
(9 = 7r/4) profiles of both the pressure field and specific charge 
distribution are presented in Fig. 12 Again, there are toroidal 
off-equatorial structures located symmetrically above and un- 
der the equatorial plane, having cylindrical distribution of the 
specific charge. 

3.2.3. Family IV off-equatorial tori 

As in the case of equatorial family IV tori, we choose 
Ki = -1. The specific charge distribution, pressure field and 
orbital velocity field are thus given by relations (|47]l-(49i. 
Further, we choose K2/\/ GM = 2.5 and 6^ = ±7r/4, giving 
the radial coordinate of tori centers = 18.75. Subsequently, 



liC/y/GM = -375. For Pq = -0.02, the meridional sections 
and off-equatorial (9 = tt/4) profiles of both the pressure 
field and specific charge distribution are presented in Fig. [13] 
There are toroidal off-equatorial structures located symmetri- 
cally above and under the equatorial plane, now having a ra- 
dial distribution of the specific charge. 

4. POLYTROPIC FLUID 

For any perfect fluid, the basic set of partial differential 
equations, describing the axisymmetric motion of charged 
fluid in spherical gravitational and dipolar magnetic fields, has 
the form of Eqs. (fTTh and ( 12 1, where = Qm(r, 9) in general. 
Further, we definethe function 



h(r,9): 



Qn 



(58) 



In the case of an incompressible fluid {g^n = const), the given 
set of equations can be written as 



dh 
dr 



= A(r,9\ 



(59) 
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Figure 10. Negatively chai'ged equatorial torus {family IV). The meridional section and equatorial profile of the pressure field (left) and the specific electric 
charge (right) are shown. 

dh 



de 



= B{r,e). 



(60) 



Now, we take a polytropic fluid with the equation of state 

P = KbI. (61) 



It is easy to show that, in this case, Eqs. 
form 



11 1 and (12 1 take the 



dh 7 - 1 
dr 7 
dh 7 - 1 



Air, 9), 
B(r, 9). 



(62) 
(63) 



equatorial and off-equatorial toroidal structures also in poly- 
tropic case. 

5. PSEUDO-NEWTONIAN MODEL 

The formalism given above can be used also for modeling 
charged toroidal structures in the vicinity of compact objects 
like neutron stars or black holes within the pseudo-Newtonian 
approach. In this case, we can use the well-known Paczyiiski- 
Wiita (PW) potential instead of the Newtonian one to describe 
the gravitational field near compact spherical body. 



Clearly, in the case of a polytropic fluid, solutions for the 
function h differ from those for an incompressible fluid just by 
the multiplicative factor (7- l)/7, where 7 is the polytropic 
index. Since in the polytropic fluid surfaces of P = const coin- 
cide with the surfaces of g,^^ = const, it is sufficient to analyze 
the function h in searching for stationary toroidal structures. 
Therefore, the discussion of incompressible perfect fluid tori 
presented in the previous section is fully relevant also for 
polytropic perfect fluid and reveals the existence of both the 



The PW potential is given by the relation ( |Paczyhski & Wi- 
|ita|1980D 

GM 

1>PW = -;3^, (64) 

where Rg = GM/c^ is the gravitational radius of the compact 
object. The PW potential is able to capture some important 
properties of the gravitational field of compact objects fol- 
lowing from General Relativity (e.g., locations of marginally 
stable and marginally bound circular orbits), which are cru- 
cial for the orbital motion of matter in their vicinit y. For 
its 'derivation' and following discussion, see ,Abramowicz| 
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P(i-, ^/4)/p, 



-0.015 




( |2009l ). 

In terms of the PW potential, the fluid dynamics is de- 
scribed by a slightly modified condition of hydrostatic equi- 
librium. In the radial direction, it takes the form 



dP 



GM 



V 

r 



sinfl 



(65) 



Figure 11. Negatively charged off-equatorial tori (family II). The meridional section and off-equatorial (9 = 7r/4) profile of the pressure field (left) and the 
specific electric charge (right) are shown. 

is also a local minimum of the pressure, corresponding to the 
cusp in the pressure distribution, i.e., to the point in which 
one of the isobars crosses itself. In dependence on the rota- 
tional law, = v^(i), this isobar can be marginally closed, 
forming the cusp, or it can be open. In the former case, little 
overflow of such isobar enables an outflow of matter from the 
torus through the cusp. 

Next, we present illustrative examples of both the equato- 
rial and off-equatorial toroidal structures constructed with the 
use of PW potential, showing the existence of cusps, but not 
presenting a detailed analysis of all the possible cases. 

In order to have some connection with previously presented 
Newtonian tori, we choose family 11 tori with Ki = -1, i.e., 
with a uniform distribution of the specific angular momentum, 
£{r, 9) = const. The corresponding specific charge distribution 
is given by relation ( [40| and the pressure field is described by 
the formula 



while in the latitudinal direction, the condition retains the 
form (|9]). 

Again, we can assume relation between the charge density 
£»e and mass density in the form ( [T0| . Since the results 
for a polytropic fluid are given as the scaled version of those 
for an incompressible fluid, we restrict our attention to the 
incompressible fluid only. 



The integrability condition satisfies the same form ( 14 1, 
leading again to the orbital velocity ^TE[ and to four possi- 
ble families of specific charge distribution ([T7[)-([20l). Analy- 
sis of necessary and sufficient conditions forme existence of 
pressure maxima, VP = and conditions ( [23] l, reveals that in 
the case of PW potential, the maxima can also be located in 
the equatorial plane (9^ = tt/2) or out of the equatorial plane 
(6*0 y 7r/2). However, apart from the local maximum, there 



GM 1 2 2 

P(r, 0) = Po + — ^ - T QmVu, + 7^ 6mK2fJ-C 

r—2R„ 2 ^3 



where the orbital velocity v^(r, 9) is given by formula (42i 



/ sin^ 61 \ 

\~) ' 
(66) 

■ formula (42 1. 
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Figure 12. Negatively ciiarged off-equatorial tori (family III). The meridional section and off-equatorial 
specific electric charge (right) are shown. 



: 7r/4) profile of the pressure field (left) and the 



In dependence on the constants K2 and C, the torus is lo- 
cated in or out of the equatorial plane. Instead of C, the lo- 
cation of the torus center can be prescribed. Note that in the 
chosen case of uniformly distributed specific angular momen- 
tum, £{r, 6) = K2. In the case of equatorial torus, we choose 
K2/V GM = 3.8 and the position of torus center = 10, which 
imply fiC/^/GM = -0.0986. In the case of off-equatorial 
torus, we choose K2/\/GM = 2.35 and the latitude of the torus 
center 9^ = 7r/4, giving = 11.16 and ^iC/VCM = -1.989. 
Constant Pq in both the cases is chosen to obtain for the crit- 
ical self-crossing isobar the zero level. The meridional sec- 
tion through the pressure field of both the equatorial and off- 



equatorial toroidal structures is shown in Fig. 14 



6. CONCLUSIONS 

In this paper, we presented a Newtonian description of 
charged perfect fluid, i.e., a perfect fluid whose particles carry 
electric charges, orbiting in spherical gravitational and dipo- 
lar magnetic fields. The equation of motion was analyzed with 
the aim of finding stationary toroidal structures — electrically 
charged thick disks, existing in prescribed external gravita- 
tional and magnetic fields. 



From an astrophysical point of view, our investigation can 
provide insight into the interplay of gravitational and elec- 
tromagnetic forces determining the vertical structure of dusty 
tori composed of electrically charged grains (dust) in which 
pressure gradients coming from mutual interactions are not 
negligible. We adopted a toy-model approach which is ob- 
viously too simple to describe real cosmic structures, such 
as those around supermassive black holes in the nuclei of 
galaxies. However, we suggest that our model demonstrates 
an interesting possibility that a small electric charge devel- 
oped on the grains affects significantly their motion through 
gravitational and pervasive (large-scale) magnetic fields. For 
the former, we assumed the approximation of central Newto- 
nian gravitating body. For the latter, we adopted an organized 
(dipolar) magnetic field with its axis oriented perpendicularly 
to the orbital plane. 

It was shown that electric charge distributed in the fluid 
changes dramatically the properties of possible toroidal struc- 
tures. Namely, apart from equatorial toroidal structures, the 
off-equatorial tori can also exist due to electromagnetic in- 
teraction of the disk with an external dipolar magnetic field. 
For the given orientation of the magnetic field (// > 0) and 
the positive rotation of the fluid {v^ > 0), the off-equatorial 
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Figure 13. Negatively charged off-equatorial tori {family IV). The meridional 
specific electric charge (right) are shown. 
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section and off-equatorial {9 = n/A) profile of the pressure field (left) and the 
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Figure 14. Negatively charged Paczyiiski-Wiita equatorial and off-equatorial tori {family II). The boundary of each tonjs is given by a marginally closed isobar 
forming the cusp which enables accretion onto the central object represented by a black circle of Schwarzschild radius. 
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tori can be only negatively charged, while the equatorial tori 
can be both positively or negatively charged. Therefore, an 
appropriate adding of both the solutions (even for the same 
specific angular momentum distribution) can give an electri- 
cally neutral configuration composed of the equatorial pos- 
itively charged torus and two negatively charged tori above 
and under the equatorial plane. 

Another interesting and qualitatively new result is the pres- 
ence of marginally closed isobars with cusp(s) in terms of 
Newtonian gravity. Note that in terms of Einstein's gravity 
(General Relativity), the existence of cusps in the structure of 
isobars is a natural phenomenon related to specific properties 
of motion in curved spacetimes, leading to the existence of 
accreting thick disks ( Abramowicz et al.]|1978 1. In the non- 
relativistic (Newtonian) case, the cusps are related to the elec- 
tromagnetic interaction between the disk with charges and 
external magnetic fieldPl In equatorial structures, the cusps 
enable accretion onto the central object (disks belonging to 
the families 111 and IV) or outflows from the torus above and 
under the disk (the family 1 disks). However, all outflows 
through the cusps are driven by the same mechanism follow- 
ing from a violation of hydrostatic equilibrium due to small 
overflow of the critical isobar forming the disk surface. 

A possible future extension of our work could be in the 
analysis of charged tori immersed into large-scale (galactic) 
magnetic fields. Such analysis could be further extended in 
order to include a cosmic repulsion. It has been shown that 
the cosmic repulsion, represented by positive cosmological 
constant, changes substantially the structure of isobars near 
the so-called static radius (Stuchlik & Hledik 1999) , where 
the cosmic repulsion is balanced by the gravitational attrac- 
tion of central bodyj^ This leads to the existence of cusps at 
the outer edge of toroidal configura tions and enables outflows 
of mat ter from such structures (Stu chlik et al.|2000t|Slany~&| 
[Stuchlik 2005 ; Kucakova et al. 201 1]l. The combined influ- 
ence of large-scale magnetic fields and cosmic repulsion can 
be well studied in the pseudo-Newtonian framework by the 
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use of cosm ological Paczyiiski-Wiita potential ( [Stuchlik & 
|KovaF2008 ), reflecting with high precision properties of fluid 
tori (Stuchlik et al. 2009) and interacting objects on galactic 
scales dS'tuchlrk & Schee|201 l ||20T2l l. 
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